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Abstract. We calculate the quantum statistical force acting on a partition wall 
that divides a one dimensional box into two halves. The two half boxes contain the 
same (fixed) number of noninteracting bosons, are kept at the same temperature, 
and admit the same boundary conditions at the outer walls; the only difference is 
the distinct boundary conditions imposed at the two sides of the partition wall. 
The net force acting on the partition wall is nonzero at zero temperature and 
remains almost constant for low temperatures. As the temperature increases, the 
force starts to decrease considerably, but after reaching a minimum it starts to 
increase, and tends to infinity with a square-root-of-temperature asympotics. This 
example demonstrates clearly that distinct boundary conditions cause remarkable 
physical effects for quantum systems. 
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1. Introduction 



Quantum systems in less than three dimensions enjoy an increasing popularity and 
importance in physics. These seemingly simple systems (rings, boxes, "dots", etc.) ex- 
hibit unexpectedly interesting properties [1] , many of which have originally been found in 
connection with quantum field theories. The recent developments in nanotechnology now 
make it possible to manufacture such quantum mechanical devices, allowing us to study 
these various features in laboratories [2-4]. With tunable control parameters — such as 
varying the magnetic Hux that is driven through a ring, or influencing the parameters 
that characterize the boundary conditions at the walls of a box — these properties be- 
come governable, giving us thus controllable quantum devices that may be useful for many 
technological applications. 

In fact, the dependence of the physical properties on the control parameters is quite 
strong. For example, the energy spectrum of a box or a ring with junctions changes 
considerably with different boundary/fitting conditions applied at the walls/ junctions [5]. 
The aim of the present paper is to demonstrate a physical consequence of this effect in 
the context of quantum statistical mechanics. Namely, we show by the example of the 
box systems how different quantum statistical behavior emerges from distinct boundary 
conditions and the corresponding distinct spectra. One of the remarkable results is that 
not only the low-temperature behavior is sensitive to the difference in spectrum but also 
the high-temperature properties depend clearly on it. 

The setting we consider is a one-dimensional quantum well/box, divided into two 
halves with an internal partition wall. Both halves contain the same number of noninter- 
acting particles. We impose the same (Dirichlet) boundary conditions at the two ends of 
the box, while the boundary conditions at the two sides of the partition wall are chosen to 
be different (Dirichlet from the left, Neumann from the right). From the different emerging 
spectra in the two half regions, we calculate the quantum statistical force/pressure acting 
on the internal wall from the left and from the right, and the corresponding net force 
acting on it, as the function of temperature. We present both numerical results (which 
are obtained by appropriate truncations of the arising infinite sums) and analytical ap- 
proximations, the latter ones aiming at understanding the low and the high temperature 
regimes. 

The number of particles on both sides, N, is arbitrary, and is not necessarily macro- 
scopically large. Our numerical results are presented for N = 100 , which is a realistic 
population number in nanoscale quantum experiments [4] . The particles are considered as 
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bosons in this paper, however, we mention that the results prove to be quahtatively (and 
partly quantitatively) similar for fermionic particles, too [6]. 

2. Low temperature regime 

The system we consider is formulated as a one-dimensional quantum well of width 21 
given by the interval [—1, 1] with a partition placed at the centre x = 0. At the end walls 
of the well the states are supposed to obey the Dirichlet boundary condition ip{±l) = 0. 
At the centre, we assume that the partition imposes distinct boundary conditions on the 
left and the right, our choice is to impose the Dirichlet one for a; = — and the Neumann 
one for x = +0: 

V'(-O) = 0, V'(+0) = 0. (2.1) 

The two half wells seperated by the partition then admit the energy levels — e^£, 
n = 1, 2, 3 . . ., with 

et = {n-l)\ e-=n^ S = '^{^f. (2.2) 

Suppose that we put N identical bosonic particles into each of the two half wells. 
The particles will then distribute among the eigenstates according to the Bose-Einstein 
statistics, 

JV = EJVn^ ^i = ^.li_, ' (2.3) 

where we have introduced h — (38 {(3 = 1/kT). Note that are determined by the 
particle number constraint (the first of (2.3)) and are dependent on the temperature. The 
forces (or pressure) acting on the partition from the right and the left are then given by 

^^ = -E^^n- (2-4) 

n 

For simplicity, in what follows we use the dimensionless force and temperature defined by 

with which the net force on the partition becomes 

A/ = /--/+, /^ = E^n<- (2-6) 

n 
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Figure 1. The net force A/(t) for N = 100, in the temperature region t < 1 , 
obtained by a numerical computation (sohd hne), and approximated by Eq. (2.8) 
(dashed hne). 

Our objective is to find how the net force A/ behaves as a function of the temperature 
variable t. 

To proceed, we first study the low temperature regime where a finite number of par- 
ticles remain in the ground state n = 1. This occurs if + bef -C 1, and in this regime 
a numerical computation exhibits interesting behaviors of the net force A/. Namely, as 
shown in Fig.l, A/ starts off by the value A/ = SN/A at t = and decreases slightly but 
basically stays there for t < 1. Above t ^ 1 the net force starts to decrease almost linearly 
until it reaches its minimum at around t ~ AT (see Fig. 2), and from there it increases to 
infinity. To see how these behaviors arise, let us consider the case of the extremely low 
temperatures t < 1, where most of the N particles are in the ground state. In this case, if 
we write in (2.6) as 

oo 

f^ = efN + J29t ^n=K-efX, (2.7) 

n=2 

then we find that the first term gives the zero temperature value /^(O) and the rests 
represent the corrections from the higher energy levels at finite temperatures. Since 
decreases exponentially fast for higher n, we may keep only the first contribution n = 2 to 
get 

A/(0 ~ ^iV+(3e-3/*-2e-2/^). (2.8) 
This gives a good approximation for t < 1 as seen in Fig.l. 
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The linear decrease of the net force observed for temperatures higher (but not much 
higher) than t=l may be understood heuristcially as foUows. First, we classify the energy 
levels into three categories, where the first is those levels for which b{e^ — ^t) <^ 1) the 
second is those for which b{e^ — ^f) ^ 1; and the third is the rest, i.e., those for which 
b{e^ — ef) is not far from 1. Let m be the level whose b{e^ — ef) is closest to 1, that is, 
6(e^ — ef) « 1. Levels with n <^m thus belong to the first category, and for these one has 
+ be^ b{e^ — ef) and hence ^^^^i-^. The corresponding given in (2.7) 

are then found to be | . Thus we have g~ — g:^ !v showing that the levels in this 

category do not contribute to the net force. On the other hand, levels with m belong 
to the second category, and for these we have A^;^ ~ ^-Ke-t-^f) but since h{e^ ~ ^f) ^ 1? 
the corresponding force is seen to be exponentially small and can be ignored. Thus 
the contribution to the net force can comes only from the intermediate levels belonging 
to the third category. As a representative, let us choose m such that ^>(e+ — e^) = 1 
holds in effect, which is equivalent to w? — m — t, for the right half well. Since we have 
Ar+ = l/(e-l), weget 



In contrast, for the level m for the left half well, we have b{e^ ~ ) = ^ t ^ ~ ^''^ m 
hence A^^ = l/(e^"'"^/"^ — 1). Thus the contribution to the force becomes 

S^=r ei+vAi - (2.10) 

Combining (2.9) and (2.10), and estimating roughly that the total number of levels in the 
third category is of the order of m, ranging from m/2 < n < 3m/2, say, we just multiply 
the number m to each of the contribution to find 

- + E iSn - 9^) = In- (2.11) 

The linear decrease of the net force is now seen in (2.11). Our argument assumes a^-\-bef < 
b(e^ — ef) which gives the upper limit of the temperature t for which the heuristic formula 
(2.11) is available. The numerical result shows that the limit is around t ~ 2N/3. 

If the temperatures are even higher t » 1 but still satisfy + bef <^ 1 remaining in 
the low temperature regime, the above approximation becomes worse and we need to resort 
to more systematic means to reproduce the numerical results. One such method valid for 
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Figure 2. The net force A/(i) for N = 100 and in the temperature region 
< i < 160, obtained by a numerical computation (solid line), and approximated 
by Eq. (2.11) (dotted line), and by Eq. (2.16) using Eq. (2.18) (dashed line). 



t 3> 1 is to consider the infinite sums as the trapezoidal approximations of corresponding 
integrals. Thus we can write 



n=l 



yjsi) 1 

2 As 



y{s)ds 



(2.12) 
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for functions y{s) vanishing at infinity, limg^oo y{s) = 0. The summation in (2.12) is taken 
at equidistant points, s^, n = 1, 2, . . ., with As = Sn+i — Sn = constant, and obviously the 
approximation is better for smaller As. Using this approximation, we obtain 



^ CO 

-T- 



Na- 



b — ! e"- 

n=l 



Na^ 1 
— +6 



±+6e± _l b 

1 a± + (s±)2 



1 \ - Q 



a± + (4) 



±^2 



b — ! e^ 

n=l 



+ is^)^ - 1 



1 



2e«^+(«f)'-l VbJsf e«*+*'-l 



-ds 



where we have introduced Cn which are equidistant on account of (2.2). Note that 

the increments As='= = Vb = are indeed small for t » 1. Since 



e^-l 



1 for small z, 
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we can approximate the second term by ^, and similarly the integral J^^ ^a±+~s^ 



st. We thus obtain 



-Na^ + 1/2- sf/Vb 



+ 637./ ^.^d.. (2.13) 



Employing the formula ^j^j ^ e ^^l + f + ff) which is handy for evaluating the integral 
approximately, we find 

r ^±^d5 ^ ^(63 - 35a±) + O (a^^) , (2.14) 

and hence 



/± ^ (^-Na^ + 1/2 - t + ^(63 - 35a±)t3/2 . (2.15) 

The net force is then 

A/ - (iVt + I^Ft^/^) K - «") + ( - V^) ^ • (2-16) 

In passing we remark that this formula (2.16) turns out to be valid even for small temper- 
atures. This can be seen, for example, by considering the t ^ limit, where ^ 
and hence from (2.3) we have ~ ln(l + 1/N) — bef . Thus the net force (2.16) is found 
to be A/ ^ iV (ef - ef) = 3iV/4 for t < 1 as seen in (2.8). 

To make use of (2.16), we need to determine the as functions of t from the total 
number costraint in (2.3). This can be done with the help of the same approximation 
method as used above. Namely, upon using (2.12) we write down the constraint condition 
as 

00 00 ^ 

^ = E = + E 

n=l n=2t^ /2J7) 

1 1/2 I r ds 

+ — ; — + ' 



where we have kept Nf" separately to achieve a better approximation (since are rapidly 
decreasing functions of n for small n). Limiting the range of integration to [0, -\/2] which 
provides the main contribution to the integral, and using -^rui ~ z ~ I valid on the range, 
we find 

1 1/2 3 1 V2 - a± - st 



a± + bef a^ + bet Vb 2 

(2.18) 



+ 



\/b\a 



±1 



where the function A is the arctan function for positive , and is the arctanh function for 
negative a"^. Unfortunately, it is not easy to solve (2.18) directly for a"^ even approximately 
to obtain an analytic formula that can be used in (2.16) to reproduce the numerical result. 
Nevertheless, the formula (2.18) has allowed us to evaluate the infinite sum to a good 
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Figure 3. The net force A/(i) for N = 100 , obtained by a numerical computation 
(solid line), and approximated for high temperatures with Eq. (3.7) (dashed line). 



accuracy, and we can solve (2.18) indirectly for a"^ by numerical means. The outcome of 
this semi-analytic analysis is in good agreement with the numerical computation as shown 
in Fig. 2. 



3. High temperature regime 

Let us now determine the high-temperature asymptotic behavior of the net force. For 
increasing temperature, we expect to decrease so [see (2.3) for n = 1] we expect 
to increase to higher positive values. Inspired by this, let us expand in := e"""*" as 



= 

n 



q e 



1 — q^e 

which is valid for any positive a^. Thus 



(3.1) 



n=l 



k=l 



n=l 



k=l 



± 1 °° 

— +2 ^ 



-kbe^ 



(3.2) 



with the constants a'^ = , a~ = 1 corresponding to the ± half wells, where we have 
extended the meaning of the notation [see (2.2)] to negative ns, too. Applying now the 
Poisson summation formula 



oo oo °^ 

E y(^)= E / dsy(s)e' 



(3.3) 



m= — oo 



we obtain 



k=l 



m——oo 



with = . Similarly, for the force f^, one can find 

oo oo I oo / 

n=l k=l m=— oo ^ 



1 - 1 e-^"" 



(3.4) 



(3.5) 



For the high-temperature asymptotic behavior (g^ ^ 0), it suffices to consider only 
the first some terms in the sums over k [both in (3.4) and (3.5)], and within each term 
to keep only the m = term in the sums over m (the m ^ terms being exponentially 

1/2 

suppressed). Now, the leading, k = 1 term in (3.4) gives that = 2N {^y' +0{b) . 
Since this leading behavior of is independent of cr^, inserting it into (3.5) gives that the 
leading, O {b~^) term of (coming from k = 1 , m = 0) is also a-independent. Hence, 
this term will drop out from the net force. Therefore, to have the first nonvanishing term 
in the net force wc need the first subleading term in g, too. Incorporating the k — 2 term 
as well for g^, we find 



2N 



1/2 



+ 2N 



- V2N 



+ 



J TT 



(3.6) 



Plugging this into (3.5) and then calculating the net force yields 



1/2 



+ o (t°) . 



(3.7) 



We can see in Fig. 3 how the net force actually reaches this square- root asymptotic behavior 
at high temperatures. 



4. Discussion 

We have found that the net force acting on the separating wall is nonzero at low 
temperatures, being practically constant for very small temperatures and starting to de- 
crease when temperature is increased. Knowing that the energy spectrum is different on 
the two halves this property is not very surprising. What is surprising, however, is that 
this decrease stops at a certain temperature and the net force starts to increase above this 
value. Furthermore, a remarkable fact is that this increase does not stop nor converges to 
some finite high-temperature limit but increases to infinity, as the square root of the tem- 
perature. From the naive expectation that such quantum effects coming from the different 
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boundary conditions should vanish at high temperatures where the classical picture would 
be available, this result seems quite unusual. However, this may be understood by the fact 
that, contrary to most quantum systems, one dimensional boxes have such energy spectra 
that the level spacing is not decreasing but increasing for higher energy levels (which is 
actually valid not only for boxes with Dirichlet and/or Neumann boundary conditions but 
for all other boxes as well [5]). In other words, quantum boxes can be distinguished by 
their high-temperature behavior, too. 

We mention that the calculation presented here could be repeated for boxes with other 
boundary conditions, too. We note however that, for most box systems, the energy values 
are determined by a transcendental equation [5], and hence a certain additional difficulty 
for carrying out the calculations, especially the analytical ones, will arise there. 
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